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Abstract. The Einstein/Maxwell equations reduce in the stationary 
and axially symmetric case to a harmonic map with prescribed singu- 
larities i/j : \ E — > Hq , where E is a subset of the axis of symmetry, 
and Mq is the complex hyperbolic plane. Motivated by this problem, we 
prove the existence and uniqueness of harmonic maps with prescribed 
singularities tp: \ E — > H, where E is a submanifold of R" of co- 
dimension > 2, and H is a classical Riemannian globally symmetric 
space of noncompact type and rank one. This result, when applied to 
the black hole problem, yields solutions which can be interpreted as 
equilibrium configurations of multiple co-axially rotating charged black 
holes held apart by singular struts. 



1. Introduction 

Let {M,g) be a four-dimensional Lorentzian manifold, and let F be a 
two-form on M. Consider the Einstein/Maxwell field equations: 



(1.1) RiCg--Rgg = 2TF 



1 

2 

(1.2) F = dA 

(1.3) d*F = 



where RiCg is the Ricci curvature tensor of the metric g, Rg its scalar curva- 
ture, and Tp is the energy-momentum-stress tensor of the electromagnetic 
field F: 

Tf{X, Y) = l- {ixF ■ iyF + ix*F ■ iY*F) 
(1-4) ^ . 

= ixF-iYF--\F\^g{X,Y), 

Here, ix denotes inner multiplication by the vector X, * is the Hodge star 

2 

operator, a ■ t denotes the inner product, and \a\ the norm of /c-forms. 
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These are given in local coordinates by: 

(1-5) *^/il.../i4-fc ~ ^ ''£l/l...i^fc/xi...^4_fc! 



(1.6) a ■ r = - * (a A *t) = ■^ct'^^-^'^t^i.../.,, 

(1-7) \a\^ = a-a = ^a^^--^^a,„„,^, 

where, according to the Einstein summation convention, repeated indices 
are summed from to 3, and e is the volume form of g. Note that tr Tp = 0, 
hence taking the trace of Equations yields Rg = 0. Thus, Equa- 

tion (jl.ip can be rewritten as 

(1.8) Ricg = 2 Tp. 

We are using rationalized units in which 4ttG = 1, where G is the gravi- 
tational constant. We will seek asymptotically flat, stationary and axially 
symmetric solutions of Equations Hl.l|) - H1.3|) . 

These equations reduce, using an idea originally due to Ernst (J], to an 
axially symmetric harmonic map (p: \ S — > H^, where S = A \ UjLi 
consists of the axis of symmetry A C with closed intervals Ij re- 
moved, and = SU{1, 2)/S{U{l) x U{2)) is the complex hyperbolic plane, 
see |121 IH]. Each interval Ij corresponds in {M,g) to a connected compo- 
nent of the event horizon. Thus will denote, throughout this paper, the 
number of black holes. The space has negative sectional curvature, and 
finite energy harmonic maps into spaces of negative sectional curvature are 
well understood, see jl7l I18j . Nevertheless, this problem is analytically in- 
teresting due to the fact that the boundary conditions for ip on T, and as 
r — > oo are singular, and hence ip has infinite energy. For the case = 1, 
there is a family of solutions known in closed form, the Kerr-Newman black 
holes, see Pj. We note that the problem has 4:N — 1 parameters: N masses 
nij, N — 1 distances dj, N angular momenta Lj, and N electric charges qj. 

Special cases of the Ernst reduction have been used to prove non-existence 
results for the Einstein equations. In early work along this line, Weyl inves- 
tigated the vacuum static case and showed that the equations reduce to a 
single linear equation. The case = 1 was known to have the Schwarzschild 
solution. Superimposing two such solutions, Weyl obtained new solutions 
which could be interpreted as equilibrium configurations of a pair of black 
holes. However, with Bach ^ in 1921, they showed that an obstruction 
arose as a conical singularity along the axis separating the two black holes. 
Having interpreted this singularity as the gravitational force, they computed 
its value and verified that the result was asymptotic to the Newtonian grav- 
itational force in the appropriate limit. The reduction was also used by 
Robinson [161, following work of Carter |2j, to prove that within the A^ = 1 
vacuum case, the Kerr solutions were unique. Mazur jl2|, and Bunting [3] 
independently, generalized this uniqueness result to the charged A^ = 1 case. 
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In j2()[ 1^ , we used the Ernst reduction to construct vacuum solutions 
which could be viewed as nonlinear generalizations of the Weyl solutions. 
Clearly, it is of interest to find out whether the obstruction found by Weyl 
in the vacuum nonrotating case persists also for all possible choice of the 
parameters, i.e., whether the force between rotating black holes is always 
positive. It has been conjectured for some time that this is so, see |15[ 
Problem 14], but the evidence is limited. Indeed, it is unlikely that the 
angle deficiencies can in general be computed exactly. We mention here a 
number of partial results, e.g., where the small angular momentum case 
is treated, [H], where non-existence is proved in the extreme regime, and 
also l^l, where Y. Li and G. Tian proved non-existence in the case when the 
solution admits an involutive symmetry. 

In the present paper, we begin the generalization of this work to the 
Einstein/Maxwell Equations. We prove the existence of a unique (4A^ — 
l)-parameter family of solutions to the reduced problem. Our main tool 
is the study of harmonic maps with prescribed singularities into classical 
globally symmetric spaces of noncompact type and rank one j23j . i.e., the 
real, complex, and quaternionic hyperbolic spaces, see 1131. However, the 
results of 23^ do not apply directly to the problem considered here. In , 
we restricted our attention to bounded domains 0, C M", and to singular 
sets S compactly contained in fi, while here O = M"" and S is unbounded. 
The necessary generalization is achieved in two steps. First, we consider 
the problem on a large ball B C M", allowing, however, the singular set E 
to meet dB. Then, we let B exhaust M". In this second step, we need the 
boundary data at infinity in R" to be given by a harmonic map which admits 
those singularities which lie along the unbounded components of S. In the 
application to black holes, this map is obtained from the Kerr-Newman 
solutions. 

The main motivation for this generalization is an attempt to come to 
a better understanding of the force between rotating black holes, i.e., the 
magnitude of the conical singularities on the bounded components of the 
axis. An important question in this respect is the dependence of this force 
on the parameters. It should be pointed out that the force can vanish in 
the case of the Einstein/Maxwell Equations. Indeed, even in the static case, 
equilibrium can be achieved with masses and charges being equal |1H I14j. 
However, these configurations are extreme, i.e., all the event horizons are 
degenerate, see Section 1231 Consequently, one could still pose the question 
of whether the force is positive for all nondegenerate Einstein/Maxwell black 
holes j^. In a future paper, we will study the regularity properties of these 
maps. This is a necessary step in order to complete the application of the 
results presented here to black holes. 

The plan of the paper is as follows. In Section [2 we review the Ernst 
Reduction as it applies to the Einstein/Maxwell Equations. We carry out 
the reduction in two steps. First in Section 12.11 we reduce the equations 
in the presence of one Killing field. Then, the second reduction, with two 
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commuting Killing fields, is carried out in Section 12.21 In the rest of Sec- 
tion |2l we examine the parameters of the problem, derive singular boundary 
conditions, and set-up a Reduced Problem for the stationary and axially 
symmetric Einstein/Maxwell Equations. In Sectional we prove our main 
existence and uniqueness result. In Sectional we address briefly the problem 
of going back from a solution of the reduced problem to a solution of the 
Einstein/Maxwell Equations. 

2. The Generalized Ernst Reduction for the 
Einstein/Maxwell Equations 

In this section, we describe the generalized Ernst reduction for the Ein- 
stein/Maxwell equations in the stationary and axially symmetric case. The 
main result in this section is Theorem^ based on which we state the Re- 
duced Problem for the stationary and axially symmetric Einstein/Maxwell 
Equations in terms of harmonic maps with prescribed singularities. The 
computations are carried out mostly in the exterior algebra formalism which 
is particularly well suited to the Maxwell Equations. 

2.1. Axial Symmetry. We first describe the connection between the ax- 
ially symmetric Einstein/Maxwell Equations and harmonic maps into H^. 
We note that the arguments given here could be applied to any spacelike 
Killing field. 

Definition 1. Let (M, g) be an oriented 4-dimensional Lorentzian manifold, 
and let F be a two-form on M. We say that (M, g, F) is SO {2) -symmetric if 
SO{2) acts effectively on (M, g) as a group of isometrics leaving F invariant. 
The axis in an 50(2)-symmetric spacetime, is the set of points fixed by the 
action of S0{2). We say that {M,g,F) is axially symmetric if it is 50(2) 
symmetric, and it has a nonempty axis. 

Let (M, g, F) be a simply connected 5'0(2)-symmetric solution of the Ein- 
stein/Maxwell Equations. Let ^ be the Killing field generator, then L^F = 0. 
Note that if {M,g) is causal, i.e., admits no closed causal curves, then ^ ei- 
ther is spacelike, or vanishes. Define the one-forms a = i^F, and P = i^*F. 
Then, we find da = —i^dF + H^F = 0. Thus, there is a function x such that 
dx = OL. Similarly d/3 = 0, hence there is a function ^ such that dil) = (3. 
Note that x ^ are determined only up to constants. Now, define the 
one-form 'y = xdip — ip dx, and observe that ^7 = 2a A /3. 

It is easy to see from ()1.5() that for any A;- form a and one- form 6, we have: 

ig*a = *{a A 6), 

where we have used the metric g to identify the tangent space TpM and its 
dual T*M. Using [£5, *] = 0, and ** = (-1)''+^ on fc-forms, it follows that 
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where 6 = *d* is the divergence operator on forms. Thus, if we define the 
twist of ^ by u; = A ^), we find 

Now for any Kilhng field ^, we have 5^ = 0, and 

(2.1) = -2trV2C = 2i^RiCg. 

Hence in view of Equation (|1.8|) . we obtain *duj = 4:i^Tp A ^. Furthermore, 
in view of (|I31), i^Tp = -i^F- {1/2) Thus, we have *duj = 4^ Aii„F. 

On the other hand, (3 = *{F A ^), hence 

*(a A /?) = -ia * /? = ^ A iaF. 

It fohows that 

(2.2) duj = 2dj, 

i.e., the one-form w — 27 is closed. We conclude that there is a function v, 
also determined up to a constant, such that 2dv = uj — 27. 

Let M' = {x e M; > 0}, and let u = -log|^| on M' . Then, we 
have i^d^ = —di^S^ + = 2e~'^^du. Thus, since 5du = —An, we have 
(5(e-2"dn) = 2e-2" \duf - e-^^An. On the other hand, i^d^ = -* (*d^ A^), 
hence using Equations ()2.1|) . (|1.8|) . and H1.4|l . we obtain: 

(5(i5de) = -hM^ + = -2(!a|^ + + • 
Furthermore, since uj = i^* d^, we find 

\ujf = *{*dC A i^d^ A e) + *(*de A A i^C) = 4e-^" [dnj^ - e"^" |(^^|2 
i.e., = 4e-2" - e^" |a;|^. We conclude that 

(2.3) Au - ^e^" Iwp - e2"(|a|2 + |/?|2) = 0. 
Now, 6uj = *{d^ A d^), and 

2 dn • w = e^" * (d^ A i^d^ A + e^" * (d^ A d^ A i^C) 
= -2dn • a; + *(d,^ A d^), 

i.e., 4dn-a; = *(dCAdO. Thus, we find that (5(e^"a;) = -4e^"du-u; + e^"5a; = 
0, which we write as 

(2.4) div(e^"a;) = 0, 

where we have put diver = —5a for any one-form a. In addition, since 
a = — * {*F A ^), we find 6a = — * d{*F A^) = d^ ■ F, and 

2du - a = -e^'^ * (d^ A % * F A ^) - e^" * (d^ A *F A 

= e^"/? -uj + di-F. 

It follows that 5(e2"a) = -2e2«du • a + e^"^ = -e^"/3 ■ i.e., 

(2.5) div(e2"Q) - e^"/3 • a; = 0. 
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Similarly, 5(e^"/?) = e^"a ■ uo, i.e., 

(2.6) div(e2"/3) + e^"a • w = 0. 

Substituting the definitions of f, x and ^ into Equations 1)2. 3^ - 1)2. 6p . it 
follows that if = (n, v, x-, "0) satisfies the following system of equations: 

(2.7) Au - 2e^" \Vv + - V'VxP - e2"([VxI^ + [VV'I^) = 

(2.8) div(e^"(Vw + - ^Vx)) = 

(2.9) div(e^"Vx) - 2e^"VV' • (Vw + xV^ - VVx) = 

(2.10) div(e2"V^) + 2e^"Vx • (Vw + xV^ - VVx) = 0. 



It is now clear that for every subset CC M', = (ti,f,x,V') is a critical 
point of the functional: 

En{v) = [\Vu\^ + e^" \Vv + xV^ - VVx|' + e2"(|Vx|' + jV^I^)} 

where d^g is the volume element of the metric g. Thus, if we choose an 
'upper half-space model' for H^, i.e., with the metric given by the line 
element: 

(2.11) ds^ = d'u2 + e^"(dz; + x#-^clx)^ + e2"((ix^ + #2), 

then the map (/?: {M',g) is a harmonic map, see Section 13 and \2'A\ . 

2.2. Stationary and Axially Symmetric Solutions. We now turn to 
the case where (M, g, F) admits an additional symmetry. 

Definition 2. Let {M,g) be an oriented and time-oriented 4-dimensional 
Lorentzian manifold and let i*" be a two-form on M. Let G be a group acting 
on M. The orbit of a point p G M is degenerate if the isotropy subgroup 
at p is nontrivial. The axis is the set of points whose orbits are degenerate. 
We say that {M,g,F) is stationary and SO {2) -symmetric if the group G = 
R X SO{2) acts effectively on {M,g) as a group of isometries leaving F 
invariant and such that the orbits of points not on the axis are timelike 
two-surfaces. We say that {M,g,F) is stationary and axially symmetric if 
it is stationary and S'0(2)-symmetric and has a nonempty axis. 

Assume that (M, g, F) is stationary and axially symmetric, and let ^ be 
the Killing field generator of the S'0(2)-symmetry normalized so that its 
orbits are closed circles of length 27r|^|. Let r be a linearly independent 
generator. Then, we have [^,t] = 0, and as before, if {M,g) is causal, then 
^ is either spacelike or vanishes, i.e., is spacelike outside the axis. We will 
now prove that if (M, g, F) is a stationary and axially symmetric solution 
of the Einstein/Maxwell Equations, then: 

(i) F and *F vanish on the orbits, i.e., F(^,r) = *F(^,r) = 0; 

(ii) the distribution of planes orthogonal to the orbits of G is integrable. 
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To see (i), note that [XL^, i^] = ij^^^j = 0, hence we have 

dira = —irda + Uri^F = 0. 

Since irCt = on the axis it follows that F{(^, r) = irO = everywhere. 
Similarly, *F(^, r) = irf3 = 0. To show (ii), it suffices by Frobenius Theorem 
to show that: 

(2.12) *(^ArAd^)=0, *(^ A r A dr) = 0. 
However, in view of (|2.2|) and (i), we have 

d * A r A d^) = dirU) = —irduj — LrU! = iiria A /?) = 0, 

where lv is the twist of ^. Thus, *(^At Ad.^) is constant, but since it vanishes 
on the axis, it is identically zero. Similarly, A r A dr) = 0. 

Let Q be an integral surface of the distribution of planes orthogonal to the 
orbits, and let h be the metric induced by g on Q. Then, the quotient space 
M/G, with its quotient metric can be identified with {Q, h). We choose the 
orientation on Q so that A r) is positively oriented. The map f. M' — > 
is invariant under G, hence reduces to a map, which we also denote by 
</5, on the quotient Q' = Q r\ M' . Indeed, let C, be an arbitrary Killing field 
generator in the Lie algebra of G. We have = 25f([C,^],0 = 0, and 

hence C,u = 0. Also, Cx = ^c*^ ~ -^(^' — 0) similarly Qijj = 0. Finally, 
7(C) =X/9(C)-V'a(C) =0, hence2C?; = ^c(^-27) = *(c??AeAC)-27(C) = 0. 

Define o" = ^ A r, then on M' we have |cjp = |r|^ - (^ • rf < 0. Let 

9 I 1 2 

p = —\a\ , then p is invariant under G, and thus reduces to a function on 
Q'. It follows that for every subset Q CC Q', the map ip: {Q', h) is a 

critical point of the functional: 

E'M = [ f|Vn|2 + e^" \Vv + xVV' - ^PVxll + e'"(|Vx|^ + pd^in, 
Jn L ^ 

where |-|^ is the norm with respect to the metric h, and d/i/j is the volume 
element of the metric h. Indeed, if O C is such a subset, and / is a 
function on M' invariant under G, then lirj^fpdfih = /c n/^/^S' '^here 
G ■ is the orbit of under G. 

Next, we show that A^p = 0. Since Agp = p^-*^ div/j(/9 V/o), it suffices to 
prove that Agp = p~^ \dp\ . To see this, note first that p^ = — <7(^,r) = 
i^irCT- Thus, we obtain: 

(2.13) 2pdp = di^ircr = i^irda = — * {*da A a). 
Therefore, we see that 

2pAgp + 2 |dpp = -6{2pdp) = *d{*da Aa) = a ■ 6da - jdcrp , 
or equivalently 

(2.14) Agp = —{a ■ 5da - |dcTp - 2 \dp\'^). 

2p 
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On the other hand, since da = /\ t — £^ f\ dr, we have, in view of Equa- 
tion (|2.12j) . that * da = *{da A ^) =0, and similarly ir * da = 0. We 
deduce, using (|2.13j) . that 

\dp\^ = *{irda A *da A i^a) = — * [da A *da A iri^a) = —p^ \daf , 

i.e., 

(2.15) \daf = -4 \dpf 

Furthermore, we claim that a ■ 6da = 0. Indeed, in view of Equation (|2.1() : 
Sda = *d{ir * d^ — i^* dr) = 6d^ A r + ^ A Sdr = 2 {i^ Ricg Ar + ^ A v RiCg), 
and consequently, 

a-6da = iri^Sda = 2 |rp RiCg(C, ^) - 4 • r) RiCg(C, r) + 2 1^1^ Ric<;(r, r). 

Introducing a = irF, and P = ir * F, we can use Equations (|1.8|) and 1)1.4(1 
to write 

Ric3(e,e) = i«i' + i/3|' 

Ricp(^, t) = a ■ a + P ■ P 
mcg{T,T) = \af + \P\^ 
from which it follows that: 

(2.16) a-6da = 2 |rp (|q| V -4 (^-r) (a-Q+/?-/3) +2 (|ap + 
In addition, we have: 

F = ^ Aa + i/:* P 
\^\'^ * F = A P - i^* a, 

and hence, 

(2.17) ICpa- (C-r)a = 

(2.18) |eP^-(e-r)/3 = -vi5*a. 
A computation similar to the one leading to (j2.15j) yields: 

I ■ ■ |2 2 I |2 t ■ ■ o\2 2 \n\2 

\lrl^*(X\ =p \a\ , \It-1^*P\ = P \P\ . 

Thus, taking the norm squared of both sides of (j2.17j) and (|2.18|) . and adding 
the results, we obtain: 

{\df + |/3|2) - 2 lel' {^■r){a-d + P-P) + [^P \r\' {\a\' + |/3p) = 0, 

which, in view of (|2.16|) . implies a ■ 6da = as claimed. Substituting this 
result back into ()2.14|) . and taking into account (|2.15() . we obtain Agp = 
p^^ as required. 

Therefore, if dp ^ 0, the function p can be used locally as a harmonic 
coordinate for the metric h on Q' . Let z be a conjugate harmonic coordinate, 
i.e., a function on Q' such that \dz^ = \dp\'^, dz ■ dp = 0, and dp A dz is 
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positively oriented. Then, in the (p, z)-coordinate system, the metric h takes 
the form: 



habdx"-dx 



.2A 



{dp^ + dz'^), 



where A = — log \dp\. Let h = e h, then /i is a flat metric on Q' , and since 
the functional E'q is conformally invariant, we have that (/? is a critical point 
of: 



|Vnj| + e^" 



pdflr 



This gives a semilinear elliptic system of partial differential equations for 
the four unknowns {u,v,x,ip)- 

We have obtained that the metric g must be locally of the form given by 
the line element: 



(2.19) 



ds' 



-p e 



w 



dtf + e^\dp^ + dz' 



where w = —€^'^{5, ■ r). All the metric coefficients can be determined from 
the map ^p. Indeed, u is obtained directly from and we will now obtain 
equations for the gradient of w and A. These quadratures are to be integrated 
after the harmonic map system has been solved. 



Define = r + w^, and observe that r/ • ^ = 0, while |ry| 



rj ■ T 



-ph'-u; 



—e p . Furthermore, we have dw = 2wdu + e ird^ = e ijjd^. Since 
= 2 ^ A dn — e^" * (^ A tj), and ^ A = a, we find i^d^ = e^" *{u; Aa), and 
hence 

(2.20) iri^ * dw 

The operator p^^i-ri^* restricted to forms tangential to Q is the Hodge 
star operator of the metric h. Since in two dimensions, the star operator 
restricted to one-forms is conformally invariant, p~^iri^* is also the Hodge 
star operator ★ of the flat metric h restricted to one-forms. Now, on one- 
forms in two dimensions *★ = —!, hence we can rewrite Equation 1)2.20(1 as: 



e'^^p-koj. 



(2.21) dw 

To derive the equations for dA, we must now use of the components of the 
Einstein/Maxwell Equations in the pz-plane: 

Ricgidp, dp) - Ricp(52, 5^) = 2 {ig^F ■ ia.F - ig^F ■ ig^F) 

Y{:\Cg{dp,dz) = 2id^F ■ ig^F. 

A straightforward coordinate computation leads to: 



A. 



u 



(2.22) 



u 



p- P 
2p 



u„ 



u:-i)+^-{xi-xi+rp-r.) 



UpU^ + -e^^'upUz + e^'^ixpXz + i^pi^z) 
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Definition 3. We collect here a few definitions related to the notion of 
asymptotic flatness. A data set {M,g,k,E,B) for the Einstein/Maxwell 
Equations consists of a 3- manifold M, a Riemannian metric ^ on M, a 
symmetric twice covariant tensor k on M, and two vectors E and B on M, 
satisfying the Constraint Equations: 

R - kj k'^ + (tr kf = l^f + \Bf , V kij -Vjtik = 2 Sijk P 5^ 

V* = 0, V* Bi = 0, 

where e and V denote respectively the volume form and the covariant de- 
rivative of g. In addition, we require that Jgg{B,n) = for every closed 
2-surface S non-homologous to zero in M, where n is the unit normal to S. 
A 3-manifold M is topologically Euclidean at infinity, if there is a compact 
set K <Z M, such that each connected component of M\i^, henceforth called 
an end, is homeomorphic to the complement of a ball in M^. The data set 
(M, g, k, E, B) is strongly asymptotically flat if M is topologically Euclidean 
at infinity, and if there is on each end a coordinate system in 
which g, k, E, B, take the asymptotic forms: 

= + v) ^ °2(r~^/^), kij = oi{r-^/^), 

E, = ^ + o,{r~'/'), B, = o,{r-^l'), 

where m and q are constants. Here, a function / is Ok{r~^) if f = o{r~^~^) 
for J = 0, ... , k, and r = A spacetime is globally hyperbolic if it ad- 

mits a Cauchy hypersurface, i.e., a hypersurface which intersects every inex- 
tendible causal curve exactly once. Given a data set (M, g, k, E, B) which is 
strongly asymptotically fiat, the Cauchy Problem for the Einstein/Maxwell 
Equations, consists of finding a globally hyperbolic spacetime (M, g), a two- 
form F, and an embedding l: M ^ M such that (M,g,F) is a solution 
of Hl.l|) - H1.3|) . g and k are respectively the first and second fundamental 
forms of L, and such that E = irE, and B = i^ * F, where r is the future 
pointing unit normal of M in M. The triple (M, g, E) is called the Cauchy 
development of the data (M ,g,k, E, B). A triple {M,g,F) will be called 
asymptotically flat if it is the Cauchy development of strongly asymptoti- 
cally fiat data. A domain of outer communications in an asymptotically fiat 
spacetime {M,g) is a maximal connected open set O C M such that from 
each point p O there are both future and past directed timelike curves to 
asymptotically flat regions of {M,g). An event horizon is the boundary of 
a domain of outer communications. 

If (M, g) is asymptotically flat and globally hyperbolic, so is any domain 
of outer communications M in (M,g), since M is then the intersection 
of future and past sets. We will restrict our attention to such a domain. 
We may then assert that (M, g) is causal. We assume that M is simply 
connected. If in addition, {M,g,E) is a solution of the Einstein/Maxwell 
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Equations which is stationary and axiahy symmetric, then there is a unique 
Kilhng field generator r, which is future directed, and such that |r| — > —1 
at spacehke infinity. Defining now ^ and p as before, it can be shown, using 
A/ip = 0, that p has no critical points, and hence can be used as a harmonic 
coordinate on all of Q' . Furthermore, the event horizon H d M can be 
characterized by the conditions: |^|^ > and p = 0. 

If \da\ vanishes identically on some component Hq of H, then we say 
that Hq is degenerate. If a component Hq of H is nondegenerate, then 
by H2.15() . we have \dp\^ = \dz\^ ^ on Hq. We will see in Section 031 
that each component Hq of -fT is a Killing horizon in the sense that the null 
generators of Hq are tangent to a Killing vector r^o which is non-zero on Hq, 
see [21 ID- Furthermore, e^'^ \dp\^ tends on Hq to a constant kq, the surface 
gravity of Hq, defined by the equation d{\T]Q\^) = —2^0% which holds on 
Hq. Thus, Hq is degenerate if and only if = and the definition of a 
degenerate horizon adopted here coincides with the standard one, see 0. If 
a component Hq of H is nondegenerate, then Hq fl Q is an interval on the 
boundary of Q] otherwise, it is a point. Indeed, the intersection of Hq with a 
spacelike hypersurface must have the topology of a 2-sphere, see We will 
assume that the event horizon H consists of N nondegenerate components, 
which we denote hy Hj, 1 < j < N. 

Let A be the z-axis in M^, and consider Q' x SO{2) = M?\A with the flat 
metric p"^ d^) ® d4> + h. Let S = A \ UjLi where Ij are the open intervals 
on the boundary of Q' corresponding to the N components Hj of the event 
horizon H. Since > on Ij, the map (/9 can be extended continuously 
across Ij. Let x (z Ij, then since Ij is of codimension 2, one can show, using 
a cut-off as in Lemma IHl that if is weakly harmonic and has finite energy 
in a sufficiently small ball centered at x. It then follows from regularity 
theory for harmonic maps jl71 ITH] that if is smooth in that ball. Thus, 
If. \ S ^ is an axially symmetric harmonic map. We have proved: 



Theorem 1. Let {M,g,F) be a solution of the Einstein/Maxwell Equa- 
tions which is stationary and axially symmetric. Assume that M is simply 
connected, that (M, g) is the domain of outer communications of an asymp- 
totically flat globally hyperbolic spacetime, and that every component of the 
event horizon in {M,g) is nondegenerate. Define ip = {u,v,x,4'), cin-d S as 
above. Then (^: \ S ^ is an axially symmetric harmonic map. In 
addition, the metric g is of the form 1)2. 19() . and the metric coefficients w 
and A satisfy Equations ()2.21() and 1)2. 22|) . 



2.3. Physical Parameters. Let S be an oriented spacelike two-sphere in 
M. We define the mass m{S), angular momentum L{S), and charge q{S), 
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contained in S by: 

m{S) = I *dT, 

OTT Js 

4^ 2 -^^ 

The mass and the angular momentum are standard Komar integrals associ- 
ated with the Killing fields ^ and r, see and the definition of the charge 
is classical. We will always perform the integration over a submanifold of 
revolution, so that the integral can be computed on the quotient Q. Let 
1 < A: < 3, let be a /c-form and let A be a fc-dimensional submanifold in 
M which is the orbit of a (/c — l)-dimensional submanifold T C Q under the 
subgroup 50(2) C G, i.e., A = S0{2) ■ T, then we find: 

(2.26) e = -2tt / i^-O, 

J A Jr 

It is well-known that Stokes' Theorem can used in ()2.23|) - (|2.25j) to relate 
the conserved quantities of two spheres S and S' which are homologous in 
M. 

Indeed, suppose that S and S' satisfy dA = S — S' for some 3-dimensional 
submanifold A, then: 

(2.27) q^S)-q{S') = -^Jj*F = 0, 

so that the charge q{S) actually depends only on the homology class of S. 
In our situation, the second homology group is generated by spacelike cross- 
sections Sj of the components Hj of the event horizon H. Hence q{S) is 
determined by the homology class of S and by the A^ parameters qj = q{Sj), 
j = 1, . . . ,N, which we may call the electric charges of the A^ black holes. 
By ^^TM . we find: 

dip. 

In view of 1)2. 27(1 . ip is constant along each component of S. Let Sj, j = 
1, . . . , A^ + 1, denote the A^ + 1 connected components of S in order of 
increasing z. We may assume that go = V'Isi = 0) and -j/'Isj = Qj-i for 
J = 2, . . . , A^ + 1. According to Equation (|1.2j) . we also have: 

which simply expresses the absence of magnetic charge. As a consequence, 
X is constant throughout S, and we may as well assume x|s = 0. 



(2.23) 
(2.24) 
(2.25) 
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In analogy, we define the angular momenta: 

(2.28) L, = Lis,) = ^ J^^ .dC = 1^^ = 11 (dv + 7)- 

However, since the electromagnetic field carries angular momentum L{S) 
is not determined by these parameters. Indeed, if dA = S — S' and if 
A = SO{2) ■ n, then: 

LiS)-L{S') = ^ I d*di = ]^ I aAp. 

i/ A. "/ £^ 

Nevertheless, since 7 vanishes on S, it is easy to see from 1)2. 28() . that v is 
constant on each component of S. We introduce the parameters Xj = 
Jj.dv, and note that Lj = (Aj + /j)/4 where Ij = fj.^- We may assume 
that Ao = = 0, and v\^^ = Xj^i ior j = 2, . . . , N + 1. 
Finally, define the A'^ masses: 



rrij 



Equation (|2.2H) implies that w is constant on each Ij. The N constants 
Wj = w\ij may be interpreted as the angular velocities of the black holes. 
Now, in view of 1)2. 13() . we have: 

dz = -kdp = 2 * 

and therefore, 

(2.29) i^*dT = -2dz + *{dihT). 
Since = 2^ A - e^" * (^ A w), we have 

(2.30) *{di^T) = -2pi.du-woJ. 
Combining Equations 1)2. 29() and 1)2. 3U() . we obtain: 

* dr = —2 dz — wlo — 2 p-k du. 
The last term vanishes on Jj, and thus we obtain: 

- / {2dz + wuj) = pj + 2wjLj, 
4 Jh 



where 2pj = Jj dz is the length of Ij in the metric h. Again, m{S) is 
not determined by these, since the electromagnetic field carries energy. Let 
dA = S - S' where A = S'0(2) • 0, then we find: 

m{S) — m{S') = - I d{wuj + 2p -k du) 

= \ I [e^'^iW? + \li?)pdpdz + 2wahl3\ . 

For future reference, we also introduce the A*" — 1 parameters rj = Jj. dz, 
j = 2,...,N. 
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We conclude this section by verifying as mentioned earlier that Hj is 
degenerate if and only if its surface gravity vanishes. If we let rjj = r + wj^, 
then rjj is a Killing field tangent to Hj, which is null but non-zero on Hj. 
Thus, Hj is a Killing horizon. Let Hj be the surface gravity of Hj, then 

on Hj, see p. 150]. However, \r]jf = — p^e^" + e~'^'^{w — WjY, hence in 
view of (|T7T|) . 

Therefore we conclude that e^" \dp'f' Hj on Hj as claimed in Section [2. ![ 
see also (21 p. 192]. 

2.4. Boundary Conditions. In order to complete the reduction, and for- 
mulate a reduced problem for the stationary and axially symmetric Ein- 
stein/Maxwell Equations, we need to derive boundary conditions for ip on 
S and as r ^ cxD in M'^ , where r = \J ^ . Since u behaves like log p near 
S and as r — > oo, these boundary conditions will be singular. They may 
be viewed as prescribed singularities for the map (/? on E and at infinity. In 
particular, ET^i{}p) the total energy of will be infinite. 

Let uq be the potential of a uniform charge distribution of S, normalized 
so that uo + log p — > as r ^ oo in M^. The map («o, 0, 0, 0) is the solution 
of the corresponding problem when \j = qj = 0, i.e., the Weyl solution. 
Define -|- 1 singular harmonic maps = (uq, Xj-i,0,qj-i): M^\S^EI^, 
j = 1, . . . , -|- 1. Note that each maps into a geodesic 7j(t) = [t, Xj,0, Qj) 
of H^. Furthermore, the distance dist{if{x),ipj{x)) remains bounded in a 
neighborhood of Sj. Indeed, if a; G Sj is not an end point of T,j, then in 
a ball B small enough about x, the function u + log p is bounded, and the 
functions v, Xi and ■)/' are even functions of p. It is straightforward now 
using 1)2.11(1 to check that dist((/3, ipj) is bounded on B. A similar argument 
can be used at the end points. Finally, in order to get a uniform estimate on 
S, one uses the asymptotic flatness of (M, g, F), see the boundary conditions 
in [2]. This motivates the following definition [23 \ 

Definition 4. Let H be a real, complex, or quaternionic hyperbolic space, 
let Q C M" be a smooth bounded domain, let S C O be a closed smooth 
submanifold of codimension at least 2, possibly with boundary, let S' C S, 
and let ip,ip' : 0, \'E — > EI be harmonic maps. We say that (p and (p' are 
asymptotic near T,' if there is a neighborhood 17' of S' such that dist((/7, cp') G 

Thus, in this terminology, ip is asymptotic to ipj near Sj. 

Now let S = Si U Eat+i be the union of the unbounded components of S. 

\ S is the domain of a harmonic map (f corresponding to a Kerr-Newman 
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solution with one black hole, such that ip is asymptotic to near Si, and 
asymptotic to (^at+i near Sat+i. The map can be written down explicitly. 
Introduce the global parameters g, /x, m, and a defined by: 

N N N-l ^ N 

j=i j=i j=i j=i 

Here q is the total charge, 2/i is the length, in the Euclidean metric of 
M^, of the interval of the z-axis from the south pole of Ii to the north 
pole of In+1, ^ and a are the mass and angular momentum per mass of 
the corresponding Kerr-Newman spacetime. Define the Schwarzschild-type 
elliptical coordinates {s, 9) on defined by: 

= {s — m + fi){s — m — fi) sin^ 9, z = {s — m) cos 9, 

where we have assumed without loss of generality that the point z = on 
A is the center of the interval I = A \ S. Then, ip = (u, V') is given by: 
(2.31) 

u = log sin 6* + - log [s^ + 0^ + 0^ {2ms - q^){s^ + cos^ 6*) " ^ sin^ 9] 

2 /IN a^fg^s — ma^ sin^ ^) cos 6* sin^ 

V = ma cos pfS — cos 0) 7^ ^ ^— h 2ma 

+ cos^ 9 

X = —qas[s + a cos 9) sm 9 

tj) = q{s^ + cos^ 9Y^{s^ + a^) cos 9 + q. 

Furthermore, the asymptotic flatness of (M, F) implies that dist(99, (^) — > 
as r ^ 00 in M^, see 

Definition 5. Assume that S C M", and let 'p^p' : \ S — H be harmonic 
maps. We say that (p and (p are asymptotic at infinity, if dist((/?, — as 
X — > 00 in J7 \ S. 

We can now formulate the reduced problem for the stationary and axially 
symmetric Einstein/Maxwell Equations. Let AN — 1 parameters 

(2.32) (//I, . . .,HN,ri, . . .,rN-i,Xi, . . . , XN,qi, ■ ■ ■ qN) 

be given, with /Uj,rj > 0. Construct the set S = A \ UjLi-^j) where Ij 
are intervals on the z-axis, of lengths 2fij, distance rj apart, and let 
S = Si U S7V+1. Define the harmonic maps pj = {uq, \j,0,qj), and the 
map : M'^ \ S ^ IH according to Equations 1)2. 31() . We associate with the 
given parameters 1)2. 32() . the data {pi, . . . , pn+i, ^) consisting of these A^+2 
singular harmonic maps. The AN — 1 physical parameters, mj, dj, Lj, and 
qj, are determined only a posteriori. Note that the parameters 1)2.32(1 are 
nonetheless geometric invariants. 

Reduced Problem . Given a set of AN — 1 parameters as in ()2.32j) , prove 
the existence of a unique axially symmetric harmonic map p: \ S — s- Mq 
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such that if is asymptotic to ipj near Sj for each j = 1,...,A^ + 1, and such 
that ip is asymptotic to (p at infinity. 

This problem will be solved in Section |21 see the Main Theorem and its 
Corollary. Clearly, of at least equal importance is the converse: given a 
solution of the Reduced Problem, is there a corresponding solution of the 
Einstein/Maxwell Equations which is stationary and axially symmetric, and 
which is a simply connected domain of outer communications of a globally 
hyperbolic and asymptotically flat spacetime with nondegenerate com- 
ponents to the event horizon? This question will be examined briefly in 
Section 01 

3. Harmonic Maps with Prescribed Singularities 

In this section, we solve the Reduced Problem for the stationary and 
axially symmetric Einstein/Maxwell Equations. We state a somewhat more 
general result, whose proof, however, requires no additional effort. 

Let Vt C ffi", and let S be a smooth closed submanifold of M" of codi- 
mension at least 2. Let H be one of the classical globally symmetric space 
of noncompact type and rank one, i.e., either real, complex, or quaternionic 
hyperbolic space of real dimension m. For our purpose, we take a harmonic 
map ip = {ip^, . . . , ip"^) : O \ E ^ H to mean a map ip £ C°°{Q \ S; H), which 
is for each CC O \ S a critical point of the energy functional: 

Eq' = / \d(pf , 
Jn' 

where {dipl^ = Yl^=i{^k^^^k'f) energy density, and (•, •) denotes the 

metric of H. This is equivalent to requiring (p to satisfy in r2\S the following 
elliptic system of nonlinear partial differential equations: 

n 
k=l 

where Tf^ are the Christoffel symbols of H. 

Definition 6. Let C M"' be a smooth domain, and let S be a smooth 
closed submanifold of O of codimension at least 2. Let 7 be a geodesic in H. 
We say that a harmonic map ip: >E[isa Ti-singular map into 7 if 

(i) <^(17\S)C7W 

(ii) (p{x) 7(+oo) as X — > S 

(iii) There is a constant (5 > such that \dip{x)\^ > 5 dist(x, S)^^ in some 
neighborhood of S. 

We will assume that we are given + 1 disjoint smooth closed connected 
submanifolds Sj C M" of codimension at least 2, possibly with boundary. 
Let J = {1 < i < + 1 : bounded}, and let S = Uj^ j be the union 
of the unbounded components of S. For p G H, and 7 a geodesic in H, let 
dist(p, 7) denote the distance from p to 7(M), i.e., inf^ dist(p, 7(t)) . 
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Definition 7. Singular Dirichlet data for the harmonic map problem with 
prescribed singularities on M" \ S into EI consists of + 2 harmonic maps 
((^1, . . . , ipN+i^'p) satisfying: 

(i) (fj : \ S — > H is S-singular into a geodesic jj of M. 

(ii) |dist(pO) ^j) — dist(po) < C, for all 1 < < + 1, where po 
is some fixed point of H. 

(iii) If : M" \ S — > H is asymptotic to ipj near for j ^ J. 

(iv) dist((^(x), 7j) — > as X — > Sj, for j ^ J outside a compact subset of 
M". 

The following result is our main existence theorem. 

Main Theorem . Let {ipi, . . . , (pN+i, '■p) be singular Dirichlet data on M" \ 
S. Then there exists a unique harmonic map ip: M" \ S ^ H, such that for 
each j = 1, . . . , + 1, ip is asymptotic to Lpj near Sj, and such that ip is 
asymptotic to ip at infinity. 

In other words, given a harmonic map (p with the correct asymptotic be- 
havior at infinity, we can modify it to a harmonic map ip with prescribed 
singular behavior near the compact components of S. It is easily seen that 
the data (931, . . . , ipN+i-, 'f) defined in Section ITU constitute singular Dirich- 
let data on \ S. Thus, as an immediate corollary to the Main Theorem, 
we obtain: 

Corollary . The Reduced Problem for the stationary and axially symmetric 
Einstein/Maxwell Equations has a unique solution. 

We note that the axial symmetry of the solution ip follows immediately 
from the uniqueness statement in the Main Theorem. 

In |21], the case EI — EI^ corresponding to the vacuum equations was 
studied, and in [22] we proved a version of this theorem for maps <p: \ S — s- 
EI where is bounded and S is compactly contained in Q. The proof of our 
Main Theorem, which combines the ideas in both these papers, is divided 
into three steps. In step one we consider the problem on a ball B C M", but 
allow S to extend to dB. Besides some minor technical points there are no 
new difficulties in this step. In step two, we use the known map ip to obtain 
pointwise a priori bounds on d{ip,ipj) and d{ip,ip) which are independent of 
the size of B. This allows us to complete the proof in step three. 

3.1. Preliminaries. The primary tools for the proof of the Main Theorem 
are the Busemann functions on Cartan-Hadamard manifolds, see [HI El- We 
note that EI is such a manifold with pinched sectional curvatures —4 < k < 
— 1. Let 7 be a geodesic in EI, then f-y{p), the Busemann function associated 
with 7 evaluated at p E EI, is the renormalized distance between p and the 
ideal point 7 (+00) £ dM. It is defined by: 

/^(p) = lim (dist(p,7(t)) -t) . 
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It is not difficult to see that the hmit is obtained uniformly for p in compact 
subsets of H. It is well-known that is analytic, convex, and has gradient 
of constant length one. The level sets, which we denote by S^{t) = {p £ 
Elj f^fip) — ^}) called horospheres. They are diffeomorphic to M*""^. 
We will also use the horoballs defined by 'B^(t) = {p G H; f-y{p) < t}, 
and the geodesic balls ]B^(p) = {q £ M; dist{p,q) < R}. Two geodesies 7 
and 7' are said to be asymptotic if dist(7(t), 7'(i)) is bounded for t > 0. 
This is clearly an equivalence relation. The boundary of HI is defined to 
be the set of equivalence classes of geodesies in H. We will denote the 
equivalence class of 7 in H by 7(+oo). Let 7 be a geodesic in H. We 
denote the reverse geodesic t 1— > 7(— t) by —7, and we also write 7(— cxd) = 
— 7(+oo). We will introduce an analytic coordinate system in H adapted 
to 7. Let u = f-y. Let vq: S_.y(0) — > M™^-*^ be an analytic coordinate 
system on S_^(0) centered at 7(0). The integral curves of V/_^ are geodesies 
parameterized by arclength. We 'drag' the coordinate system vq along these 
integral curves. More precisely, let $t be the analytic flow generated by 
this vector field, then <I>_j maps ^--yit) to S_^(0), hence vt = voo^_t is an 
analytic coordinate system on S_^(t). Define v: M. ^ W^~^ by v\g_ = vt, 
and let $ = (u, t;) : EI — > R"*, then $ is an analytic coordinate system on H. 
In these coordinates, the metric of EI reads: 

ds^ = du^ + Qp{dv), 

where for each p E EI, Qp is a positive quadratic form on M™^"*^. We note 
that the convexity of implies that for each non-zero X G M™~^, and 
each fixed v G M™~^, Q^-i(u,v){^) is a positive increasing function of u. 

The following lemmas were proved in j23j . They are quoted here without 
proofs. Except for LemmaEl they depend only on the bound —4 < k < —1 
for the sectional curvatures of M. 

Lemma 1. For any {u,v) G and any X G M™-\ there holds: 

d 

(3.1) 2Q$-i(„,,)(X) < — {Q<i>-^(^u,v){X)) < 

Lemma 2. Let j be a geodesic in EI. Then for any to G M, and any T >0, 
we have 

'B^i-to + T) n B_^(to + T) c 1^(7(^0)), 
where R = T + log 2. 

Lemma 3. Let 7 and 7' he geodesies in H, such that j{—oo) = 7'(— 00) and 
7(-|-oo) 7^ j'{+oo). Then for some d G M, we have: 

(/-7 - fr) °7'(i) = d 

hm (/_^ + /^)o7'(t) =0. 
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Lemma 4. Let ^ he a geodesic in H. Then, for any to € M, and any 
T > (1/2) log 2, we have 

(V^(p), V/_^(p)) > 0, Vp G S^(-to + T) \ S_^(to + T) 

where (•, •) denotes the metric on H. 



Lemma 5. Let ^ he a geodesic in H. Then there is an analytic coordinate 
system ^ = {u,v) : M ^ M™" with u = such that in these coordinates, 
the metric ofM is given hy: 

ds"^ = du^ + Qp{dv), 

and satisfies the following conditions: 

(i) Let R > 0, to M, and let 7' be a geodesic in H with ^'(—00) = 
^{—00). Then, there exists a constant c > 1 such that for all t > to, 
and all p E ]B/j(7'(t)), there holds: 

\ Qj'it)iX) < Qp{X) < cQy(t)(X), VX G R"'-\ 

(ii) For all t,t' G M, the set S_^(t) fl 'B^'(t') is star-shaped in these 
coordinates with respect to its 'center', the point where 7' intersects 
S_^(t). 



The following lemma is a slight generalization of Lemma 8 in |23j . 

Lemma 6. Let B = B^i C M" he the ball of radius R centered at the origin, 
and let T, be a smooth closed submanifold ofW^ of codimension at least 2. 
Suppose that u G C°°{B \ E) satisfies Au > 0, and 0<u<linB\'S. 
Then, for any R' < R, there holds 

(3.2) / \Vu\'^ < C, 



where C is a constant that depends only on R, R' , and n. Furthermore, we 
have: 

(3.3) sup u < sup u. 

B\T, dB\T. 

Proof. If X S Cq'^(-B \ E) with < x ^ then integrating by parts the 
inequality 'x^uAu > 0, and using < u < 1, we obtain: 

(3.4) / x'|Vnp<4 / \Vxf. 

JB J B 

Define the cut-off function Xe by: 



Xe 



2 — log r / log e \i <r < e 

if r < £2 

1 if r > e, 
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where r{x) = dist(x,S), and let x' S C^{B) be another cut-off satisfying 
< x' ^ 1 and x' = 1 on 3^/. Now, take x = XeX' iii (EH])) and then let 
e ^ to get: 

/ (XO' iV^xp < 4 / \Vxf, 

JB J B 

see |231 Lemma 8], The estimate (|3.'2p follows immediately. Now, the same 
argument shows that Au > weakly throughout 5, hence, by the maximum 
principle, we obtain (|3.3() over Bp-^i instead of B. Finally, let R' ^ R to 
get (I2I31). □ 

3.2. Step One. In this step, we restrict our attention to a ball B C M" 
large enough so that (iv) in Definition [3 holds near dB. We prove the 
existence of a solution ^p: i?\S^EIina manner similar to the proof of |23| 
Theorem 1]. We will assume all the hypotheses of the Main Theorem. Some 
of the elements in the proof of Proposition ^ will be used also in the next 
steps. 

Proposition 1. Let B C M" he a large enough hall. Then there is a unique 
harmonic map ip: B \T, ^ M such that ip is asymptotic to (pj near T,j for 
each j = 1, . . . , N + 1, and p = (p on dB \ S. 

Proof. The proof follows closely that of Proposition 2 in |23j . We begin with 
the uniqueness. Let ^p and ^p' be two such maps, then clearly dist((/3, p') G 
L°°{B \ S), and Lp = (p' on dB \ S. Consider the function u = dist((/9, v'')^. 
We have that Au > on i3 \ S, see |19j . u is bounded on B \Ti, and 
u\qb\£ = 0. Thus, by Lemma El it follows that u = and hence ip = ip' on 
B\E. 

For the existence proof, we set-up a variational approach. Pick a smooth 
open cover of M" which separates the S,'s, i.e., a collection {^ j^lV of 
smooth open sets such that {jj^'j = C O^-, and n ^l'-, = 

for j 7^ j' . It can easily be arranged that if j G J then Qj CC B. Let 
Qj = ^I'jDB. Let {xj}f=i^, be a partition of unity subordinate to the cover 
{^j}f=i^ of B such that Xj = 1 near T,j. 

We may assume, without loss of generality, that all the geodesies 7^ have 
the same initial point on dM. Indeed, if this is not the case, we can write 
(pj = 'jj o Uj, where Uj is a harmonic function on B\Y,. Since IHI has strictly 
negative curvature, there are geodesies 7^ with a common initial point poo = 
7j(— 00) G dM such that 7j and 7^ are asymptotic. Let ip'j = -y'joUj, then 
it is easy to see that pj and p>'j are asymptotic near so that we may 
replace ipj by ip'j. We may also assume that the t/j's are all equal. Indeed, 
they differ at most by constants, since, by (ii) of Definition [3 Uj — Uji is a 
bounded harmonic function on So we may replace p>'j by p" = 7^- oUq, 
where Uj = uq + Cj for some -|- 1 constants Cj . Again ip" is asymptotic to 
ipj near Sj. Note that the other conditions in Definition |7| are not affected 
by these changes, and that the singular Dirichlet data (99", . . . , p'jsf^i, p) is 
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equivalent to (991, . . . , ipN+i,'f)- We will assume these changes have already 
been carried out. 

Let $ = {u,v) be the coordinate system on EI with u = f-^^ given in 
Lemmaini We will identify any map (p: B\Ti — > H with its parameterization 
$ o = {u,v) whenever no confusion can arise. We have (pj = (uq, wj) where 
Wj G M™^-*^ are constants. We write uq = X^jS^ where Uj is a harmonic 
function on which is singular only on T,j. Without loss of generality, we 
may assume that for j € J, Uj{x) = on dB^ and in particular Uj > in 

For any domain 0, define the norms: 

1/2 

(3.5) ||m 



1/2 

(3.6) 



^{b!' + Q^,(v^)}) , 



where Qipj{yv)\x = Y^^=iQ^j(x)iy kv{x)) for x E \ Sj. We will use the 
Hilbert spaces: 

Hiiyi) = {u-.Q.-* M; li-ull^ < 00} 
i/f^(17;M— 1) ={v:n^ W^'^; \\v\\^^.^^ < cx)}, 

and the subspaces Hifi{il,) and H'^q{Q;M.'^^^), defined to be the closure in 
Hi{n) of C^{n), and the closure in Hf'{n;R"'-^) of C^{n \ Sj;M™-i), 
respectively. 

We write (p = {u,v), and introduce the function vq obtained by truncating 
V near the bounded components of S and setting it to its prescribed values 
Wj € W^-^: 

Note that vq = Wj near for j G J, hence clearly vq — wj G Hf^ {^j;M."^^^). 
We also decompose uq = Uq + Uq into Uq = '^j-> the contribution from 

the bounded components of S, and C/q = ^jf^jUj, the contribution from 
the unbounded components. Clearly, AC/q = ^^^o = on i? \ S and Uq = Q 
on dB. Since (p is a. harmonic map, u is subharmonic on M" \ S, hence 
A(n — Uq) > on \ S. Furthermore, \u — Uq\ is bounded on the ball 
Brj^i. Indeed, on the one hand, for j ^ J, we have [{i — C/qI < — uq\ + Uq < 
dist(<^, (/9j) + Cj on Q.'- n Br+i for some constant Cj. On the other hand, 

for j G J, we have that both u and Uq are bounded on Hence it follows 
from Lemma El that u — Uq € Hi{B). In addition, u satisfies the equation 
A-S = {d / du){Q(p{Vv)) , and Lemma^implies that: 

Q^{Vv) < Ia{u-Uo), 
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in -B_R+i \ S. Multiplying by a cut-off x ^ i^i the proof of Lemma El and 
integrating over B, we conclude, after taking the appropriate limit, that 

(3.7) / Q^{Vv) < oo. 

Jb 

Define "K to be the space of maps (p = {u, v) : B \ ^ M. satisfying: 
\ - u - Uo e Hi^oiB); 

For maps if (z 'K define: 

F{ip) = {|V(n - no)|' + Q^{Vv)] . 

For R > 0, define Jiji to be the space of maps ip G 'K such that dist((/3, ipj) < 
R for a.e. x E Qj\T,j for j = 1,...,A^ + 1. We will need the following lemma 
which is a direct consequence of Lemma |SJ 

Lemma 7. Let R > 0, then there is a c> 1 such that for all if E there 
holds for each j = 0,... ,N and for all X G M™~^- 

^Q,p,ix){^) < Qipix){X) < cQ^^,(a,)(X), for a.e. x e 

We already know that QifijC^vo) < oo for j € J, and it follows from 

Lemma |7| and ()3.7() that the same holds also for j J. Indeed, since (f is 
asymptotic to ipj near Sj for j ^ J, there is a constant c > 1 such that 

/ Q^,{yv)<c Q^{Vv)<^. 

Since Xj is smooth and Vxj = near S, we obtain that /^^ Qipjiyvo) < oo 
as claimed. Now, let ip = (n, v) € then there is a constant c > 1 such 
that: 




Since also u — uo = {u — u — Uq) + {u — Uq) G Hi{B), we conclude that F is 
finite on 3i. 

It is straightforward to check that a minimizer 99 € IK of F is a har- 
monic map on B\T,, see j23| . Hence, by the regularity theory for harmonic 
maps ^1 E], we have ip G C°°{B \ E;EI). Furthermore, ip is asymptotic 
to ipj near for j = 1, . . . , + 1 by construction, and = (p on dB \ S. 
Thus, to prove Proposition it suffices to show that F admits a minimizer 
in "K. 

It can be shown, exactly as in |23[ Proposition 1], that F admits a mini- 
mizer in "Kr. The main steps are as follows. Consider a minimizing sequence 
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= {uk,Vk) in 'Kr. Then for some subsequence k', Uk' — uq converges 
weakly in Hi(B) and pointwise a.e. in B. Furthermore, using the bound 
dist{ipk, fj) < i? in Qj, it is shown that, perhaps along a further subse- 
quence, Vk' —vq also converges weakly in each H'^q{B,W^~^), and pointwise 
a.e. in B. Let (p = {u,v) € JC/j be the weak limit. Then clearly we have 

(3.9) / |V('u-tio)|^ < liminf / \V{uk'-uo)f. 

Jb Jb 

Furthermore, is equivalent to Q^p^ on Qj, hence we have 



(3.10) 



where 



B 



< liminf 



1/2 



Xk 



if / 
otherwise. 



However Xk' 
that 



1 pointwise a.e. in B, and Xk' is bounded. Thus, we conclude 



lim 



Xk'QA^v) = / Q^{Vv), 

B JB 



and therefore 1)3. 10(1 implies: 

j Qv(.Vv)< Q Q^iVv)^ ' liminf (^y" Q^^.i^Vk' 



1/2 



Combining this with (|3.9|) we obtain 

-^(y) ^ liminfF((^fc/ 



infF, 



and it follows that 93 is a minimizer in Kr. 

Thus, it remains to prove that for some i? > large enough, inf^F = 
inf F . This is the content of the next lemma. 

Lemma 8. There is a constant R > such that for every e > 0, and every 
ip g!K, there is (p' G !Kr such that F{(f') < F{(p) + e. 

Proof of Lemma\^ Let !K* be the space of maps ip = {u,v) G "K such that 
V = Wj in a neighborhood of Sj for j £ J, ip = (p outside some compact set 
K C B, and v = Wj in a neighborhood of n K for j J. Lemma |H1 will 
follow immediately from: (i) for each 93 G IK and each e > there is ^p' G "K* 
such that F{ip') < F{ip) + e; and (ii) for each (p G "K* there is ip' G "Kr such 
that F{(p') < F{p). The proof of (i), a standard approximation argument, 
is practically unchanged from [231 Lemma 11]. We turn to the proof of 
(ii). The bound dist{(p, ipj) < R will be achieved consecutively on each 
beginning with all j J. The case j G J is simpler. Assume without loss of 
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generality that 1 J. Introduce a 'dual' coordinate system $ = (u, v) with 
u = f-y^ as given by Lemma El Write 

^2 1-2 



ds = du + Qp{dv) 



for the metric of H in these coordinates. Set uq = — ni + Ylf=2^ '^j- V ^ 
"K* and write ^ oLp = [u,v). Then, since \ Vu^ + Qip{Vv) = \Vu\^ + Q^{Vv), 
we find that in B\T,: 

\V{u - uo)\^ + Q^iVv) = [V(n - uo)\^ + Q^{Vv) 

N+l N+1 

- 2 V(ii + u) ■ Vni - 2 ^ V(u - • Vuj + 4 ^ Vuj ■ Vui. 

i=2 i=2 

The idea is now to integrate this identity over B, use the right hand side to 
truncate u — uq, and the left hand side to truncate u — uq. The convexity 
of the Busemann functions u and u on H, expressed in the monotonicity of 
Q and Q, ensures that these truncations do not increase F. This can be 
viewed as a weak form of a variational maximum principle. 

Some care, however, needs to be taken because of the singularities on S. 
Integrating first over B\T,^ , where = {x G B; dist(x, S) < e}, we obtain: 



V{u - uo)f + Q^iVv)} = f ||V(n - uo)f + Q^{Vv)}+I,{^), 



where: 

(3.11) h{^)= I {-2div((n + u) Vui) 

N+l N+l 

- 2 ^ div((n - u) Wuj) + 4 ^ \7uj ■ Vui 
i=2 j=2 

A tedious but straightforward verification shows that Is{f) — > C as e ^ 0, 
where C is a constant independent of 93 G "K* , yielding an integral identity: 

(3.12) F{ip) = J^{\^iu - uo)\^ + 0^(V?;)} + C. 

Note that if e > is small enough, we can decompose d{B \ T,^) into a 
disjoint union Uj^i^ ^^i^j) U ("^-^ \ ^^)- Then, for example, the first term 
in (|3.11|) can be integrated by parts: 
(3.13) 

-2 div((n + n)Vni) = -2 I + /_ + E L 1 
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For the first integral on the right hand side of (jS^Sl, we have, as e — > 0: 

3.14 / (u + u)—^^ n + n — 

JaB\s^ on Job on 

which clearly is independent of (/? G !K*. The integral on the right hand side 
of ()3.14|) is finite since u + u is bounded in J7i \ Si. Next, since ip € IK*, 
there is a compact set K C B such that ip = (f outside K and v = wi in a 
neighborhood of Si n K. We find: 

dui f / , -X dui 



9Sf dn J9Y.i\K dn 

since for x G SSf n K and e > small enough, ip{x) lies along the geodesic 
7i where u + u = 0. However for x € c^Sf \ we have ip{x) = (p{x) which 
implies [u{x) +u{x)) < dist((^(x), 71) ^ as e — > 0, by (iv) of Definitional 
Thus, we obtain 

Jdi:i\K on 

Finally, in view of Lemma|31 we have that, for 2 < j < + 1, lim£_to(u + 
u)\ds = dj, for some constants dj. Hence the last sum in 1)3. also tends to 
zero since ui is regular in i}j for j = 2, . . . , N + 1. The other terms in 
are handled similarly. See 2.3^ Proof of Lemma 13] for more details. 
Now, truncate u — uq above at 

Ti = sup(n - uq) + 1, 

as 

i.e., define a map ip' = [u' ,v') E IK by ^099' = {u' ,v) where 

u' — uq = min{n — uq, Ti}. 
Clearly the map 93' E K and satisfies: 

(3.15) ip'{x) e'B^^{uQ{x)+fi), VxGS\S, 

and also, in view of 1)3. 12() . F{(p') < F{ip). Let ci = sup^^^ Yl!j=2 then 
we obtain from (|3.15j) : 

(3.16) 99'(x) G lB^,(-ni(x) +ci + ri), VxGOi\Si. 
Next truncate u' — uq above at 



Ti = max < sup(n — uq) + 1, Ti, (1/2) log 2 
I dB 

to get a map ip" = {u" ,v') G % satisfying: 

(^"(x) G S_^,(uo(2;) +Ti), VxG-B\S 
and F{ip") < F{ip'). As before it follows that: 
(3.17) (/?"(x) G S_^,(ui(x) + ci +ri), VxGJ1i\Si. 
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From Lemma m we obtain, as in j2,S| Lemma 11], that the bound (|3.16|) still 
holds for 99": 

(3.18) ip"{x) e'B^,{-ui{x) + ci+Ti), VxGJ1i\Si. 
Hence, combining (|3.17|) and (|3.18|) with Lemma [3 we conclude: 

(3.19) /(x) gBr,(v9i(x)), Vxef^iVSi, 

where i?i = ci + Ti + log 2. 

Consider now the map (p"\^jN+i^ , and notice that from (|3.19j) we can 

obtain a pointwise estimate for ip" on 

^(Uf=V%)- Thus we can proceed by 
induction to obtain a map ip'" G "K satisfying F{ip"') < F{(p), and for each 
j = l,...,N+l: 

(p"'{x) G Mr. {fjix)) , Vx G Qj \ T,j, 

for some constants Rj depending only on the data {ipi, . . . , (pN+i, 'f), Setting 
R = maxjRj, we have obtained p'" G 'Kji with F{ip"') < F{ip). This 
completes the proof of Lemma |H1 and of Proposition^ □ 

Remark. It is important to note that the a priori bounds dist((y9, (/?j) < i? in 
\ Tij given by Lemma |H1 depend on the radius of -B, hence cannot be used 
to obtain a solution on M" \ S. This is remedied in the next step where we 
obtain a priori bounds independent of the size of B. 

3.3. Step Two. In this step, we furnish the main ingredient in the proof 
of the Main Theorem. We establish uniform pointwise a priori bounds for 
d\sl{'p,(p) near 00, and for d\sl{'p, (pj) near Sj, where (/? : i? \ S — > EI is the 
solution given by Proposition^ As in Proposition^ B C M" is a sufficiently 
large ball, but we now use a slightly different open cover of M". For j ^ J ^ (p 
is asymptotic to (pj near Sj, hence we can choose Q'^ to be a neighborhood of 
Sj such that dist((^, ipj) is bounded on n'-\T,j. However, to get an open cover 

of M"", we add another open set il' which we may choose so that Sj Pi fi' = 
for all 1 < j < iV + 1. We set % = n B, and 9. = Vt'r\B. We also change 
the normalization of the harmonic functions Uj for j G J, so that Uj (x) — > 
as r = |x| ^ 00 in M". 

Proposition 2. There is a constant R> independent of B such that if ip 
is the solution given by Proposition^ then 

(3.20) dist{ip,ip) < R, inn, 

(3.21) dist{ip,ipj) < R, inflj\T.j,forj = l,...,N + l. 
Proof. Define on B\T, the function: 

1/ = \/l + dist(v9,(^)2 - Uq, 

where recall that Uq = J2jeJ ^i) defined on pagel^ is the contribution to uq 
from the bounded components of S. Note that since the function dist(-, •)^ 
is convex on EI x H, we have Au > on \ S, see |19j . We claim that i' is 
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bounded, and therefore Lemma IHl apphes. To see this, fix first j J, and 
note that (p is asymptotic to ipj, which is asymptotic to <p near T,j. Hence 
dist(v3, (^) is bounded in Qj. Also, Uq is bounded in Qj, and hence is 
bounded there. Similarly, u is bounded on Q. Now, let j € J, and observe 
that dist(7j(0), (^) is bounded on Qj. Thus, since dist (93^,7^(0)) = uj, we 
see that 

u <1 + dist{ip,ipj) + dist(99j,7j(0)) + dist(7j(0), (^) - Uq 
= 1 + dist(v9,93j) +dist(7j(0),(^) -J^j'eJUf 

Since is asymptotic to (pj near Sj, and Uji is bounded in \Sj for j' 7^ j, 
we obtain < Cj in \ Sj. Hence, we conclude ly < C in B \ T, . Similarly, 
we obtain 

u > — dist(7j(0), (^) — dist{ip, ifj) — J2j'eJ'^j' — 

Applying Lemma IHl and using the fact that ip = (f and Uq = on dB \ S, 
we deduce that: 

(3.22) u < sup z^ < 1 
It follows that 

(3.23) dist{(p,<f) < ^/{l + Uof^. 

Now there is a T such that Uq < T on Uj^j ^'j U Thus, we obtain 

immediately dist((/9, (p) < R on Q with R = 1 + T. Furthermore, for j ^ J 
there are constants Tj > such that dist{(p,pj) < Tj on fi^- \ Sj. Thus, 
from ()3.23() . we obtain the pointwise a priori estimate: 

dist((/3, ipj) < Rj, 

in ^}j\T,j for j ^ J, where Rj = 1+T+Tj is clearly independent of the radius 
of B. Using these estimates, we can now bound dist((/9, pj) on dO,j for j G J. 
Therefore, for j € J, the same argument, using the bounded subharmonic 
functions I'j = y/l + dist((^, pj^ on \ Sj, yields dist{ip,pj) < Rj, for 
some Rj independent of the radius of B. Setting R = max^ Rj, Proposition[21 
follows. □ 

3.4. Proof of the Main Theorem. The proof of the uniqueness statement 
is almost the same as in Proposition ^ If p and p' are two solutions, 
then they are asymptotic near each T,j, and asymptotic at infinity, hence 
dist(99, ip')'^ is a bounded subharmonic function on M" \ S. From Lemma 
it follows that for any ball B C M" of radius R centered at the origin: 

sup dist((/?, (/?') < sup dist((/?, (/?'), 

B\E dB\T. 

Since the right hand side tends to zero as the radius R of B tends to 00, it 
follows that dist((/9, ip') = 0, and hence ip = ip'. 



28 



GILBERT WEINSTEIN 



To prove the existence of a solution, we choose a sequence of radii Rk — 
oo, let Bk be the ball of radius Rk centered at the origin, and denote by 
'Pk = {uk,Vk) the solutions given by Proposition ^ on Bk- We will use 
the uniform pointwise a priori bounds given in Proposition |2l to prove the 
convergence of ipk to a solution on M" \ S. 

We use the open cover introduced in section fS.'Sl For any ball 

B, define the space !K'{B) as the space of maps ip: B \ 'E ^ M satisfying 

'u-uo£ Hi{B); 

disti^,ip,)eL^{n'^\^j), yj = i,...,N + i, 

and for R > denote by the space of maps (p £ 'K'{B) satisfying: 

dist((/9, (fi) < R, in Q, 

dist(v9, ipj) < R, in \ Sj, for j = 1, . . . , A'' + 1. 

Now, fix B = Bkg, let B' = Bpi^^i, and consider the sequence of maps 
(pk\B for k > kt). By Proposition |21 there is i? > such that (pk € 'K'^{B). 
It follows that on UjS^ \ '^^ have \uk — uo\ < it!. Similarly, on Q, 
we have \uk — uo\ < \uk — u\ + \u — uq\ < R' where R' = R + sup^ \u — uq\. 
Thus, we conclude that \uk — uq\ < C on i? \ S, where C = max{i?, i?'}. 
Furthermore, we have A('Ufc — uq) > on S \ S. Now the argument in the 
proof of Lemma El shows that there is a constant C independent of k such 
that 

/ |V(nfe-no)|' <C', 
Jb' 

and similarly, using Q^^{Vvk) < (l/2)A('Ufc — no), we deduce, as in the 
argument leading to (|3.7j) . that 

/ Q^d'^vk)<C". 

JB' 

Thus, we have a uniform bound F[(pk) < C" . It is straightforward now 
to show that there is a subsequence, without loss of generality pk, which 
converges weakly in 'K'{B') and pointwise a.e. in B', see |28| Proof of Propo- 
sition 1]. 

Repeating this argument for each ko, and using a diagonal sequence, it is 
clear that we can choose a subsequence, without loss of generality pk again, 
which converges pointwise a.e. in M" \ S, and which for each ball B C M"", 
converges weakly in lH'(i?) to a map p> £ !K'j^{B). For each open set O CC 
M" \ S, pk\o is a family of smooth harmonic maps with uniformly bounded 
energy into H, which maps into a fixed compact subset of H. Hence, using 
standard harmonic map theory, one can obtain uniform a priori bounds in 
C^'"(0) for ifk, and hence we deduce that a subsequence converges uniformly 
in O together with two of its derivatives. We conclude that ip is a harmonic 
map. It remains to see that ip is asymptotic to p at infinity. However, from 
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the estimate (|l-i.2:-{|) on v which holds for each (^fc, we deduce that the same 
holds of Lp: 

dist((^, (p) < ^Uq{2 + Uq) 

Since Uq{x) — > as x ^ cxo, this shows that d\si{p^ip) — > as x — > oo in 
B\Ti. This completes the proof of the Main Theorem. 

4. N- Black Hole Solutions of the Einstein/Maxwell 

Equations 

In this section, starting from a solution (/9 : \ S ^ of the Reduced 
Problem, we construct a solution {M,g,F) of the Einstein/Maxwell Equa- 
tions. 

4.1. The Spacetime Metric g. Let (/?: M'^ \ S ^ be a solution of the 
Reduced Problem. Then ip = (u, v, x, V') satisfies on ]R^\S the system (|2.7p - 
1)2.10(1 of partial differential equations, where however now, the Laplacian, 
the divergence, the norms, and the inner-product are all with respect to the 
Euclidean metric of M^. 

Introduce cylindrical coordinates (p, (j), z) in M^, and denote by ^ and Q the 
vector fields d/d(j) and d/dz respectively. Let * be the Hodge star operator 
of the Euclidean metric in and * = p~^i^* the Hodge star operator of the 
flat metric ds'^ = dp^ + dz"^ in the pz-plane. Setting uj = 2 {dv + xdip dx) 
it follows from (|2.<S|) that the two- form e^" *a; is closed on \ S, and hence, 
since w is invariant under ^, the one-form e^"i^ * u; is closed on \ E. Since 

\ S is simply connected, there is a function w, defined up to an additive 
constant such that 

(4.1) dw; = e^"ig *w = pe^"*u;. 



see (j22TI). 

Now let T^p be the stress tensor of the map p: 



T^{X, Y) = {dp{X), dp{Y)) - i \dp\' X ■ Y, 



for X,Y E M^, where (•,•) denotes the inner product in H^. Clearly is 
symmetric, and since p is a, harmonic map is also divergence free. Thus, 
since C is a Killing field in R^, i^T,^ is a divergence free vector field on R^ \$]. 
As before, it follows that * icT^p is a closed one-form on R^ \ S, hence there 
is a function A, defined up to an additive constant such that 

d{X — u) = —2i^ * icTip. 



Note that 

-2zg*i^r^ = p 



^l-ul + -/^{u^l - u:l) + e'-ixl -xl + i'l- ^l) ) dp 
+ 2 iupU:, + ^e^^'wp e^'^iXpXz + ^p^^) j dz 
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and compare with Equations (|2.22|) . It is clear that w and A defined in this 
way are axially symmetric. Furthermore, we deduce from (|4.1j) that e^"*a; = 
d(j)/\dw. Thus, Equation (|2.1flj) imphes that the two-form e'^^*dip+w dx/\d(l) 
is closed. It follows that the one-form e^^i^ * dip — wdx^s closed, and hence 
there is a function x such that 

dx = e^"z^ * dip — w dx- 

Let A C be the z-axis and let g be the metric on M' = M x (E^ \ A) 
given by the line element: 

(4.2) ds^ = -p'^e^'^dt^ + e-^''{d(l)-wdt)^ + e'^^{dp'^ + dz'^). 
Define the one-form A on M' by: 

(4.3) A = -{xd<l) + xdt), 
and let 

F = dA = dt Adx + dc/) A dx- 
We have that i^F = dx, and irF = dx, where r = d/dt. Clearly ^ and r are 
commuting Killing fields of (Af g) which generate an abelian two-parameter 
group G = M X SO (2) of isometrics leaving F invariant with timelike two- 
dimensional orbits, and hence (M', g, F) is stationary and axially symmetric. 
Furthermore, it is not difficult to check that uj is the twist of ^, and that 
* F = dip where here * denotes the Hodge star operator of the metric 
g. It is then quite straightforward to verify that {M',g,F) satisfies the 
Einstein/Maxwell equations. Finally, since the 4A^ — 1 parameters pj, rj, Xj, 
and Qj are geometric invariants, these solutions form a (4iV — l)-parameter 
family of stationary and axially symmetric solutions of the Einstein/Maxwell 
Equations. 

4.2. Conclusion. We have shown the existence a (4A'^— l)-parameter family 
of solutions of the Einstein/Maxwell Equations which are stationary and 
axially symmetric, and which should be interpreted as N co-axially rotating 
charged black holes in equilibrium possibly held apart by singular struts. In 
order to complete this interpretation, it is necessary to show that the metric 
g and the field F can be extended smoothly across the axis of symmetry S. 
This requires first the smoothness of the metric and field components across 
S, which would follow from the smoothness of e", v, x a^id ip on 9S, 
where dT, denotes the set of 2N endpoints of S; less regularity is expected on 
dTj, see |21| . For the target H = corresponding to the Einstein/ Vacuum 
Equations this was established in [20]. Regularity of harmonic maps with 
prescribed singularities into was studied in further generality in |lf)j . 

Even after smoothness is established, there is still the possibility of a 
conical singularity on the bounded components of S. As in the vacuum 
case, this is to be interpreted as a singular strut holding the black holes 
apart, and the angle deficiency can be related to the force between these. 
Finally, to prove that g is asymptotically flat, the decay estimate obtained 
at infinity must be sharpened, see [2T] . 
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31 



References 

[1] R. Bach and H. Weyl, Neue Losungen der Einsteinschen Gravitationsgleichun- 

gen, Mathematische Zeitschrift 13 (1921), 132-145. 
[2] B. Carter, Black Hole Equilibrium States, in Black Holes, edited by C. DeWitt 

and B. S. DeWitt, Gordon and Breach Science Publishers, New York, 1973. 

[3] B . Carter, Bunting Identity and Mazur Identity for Nonlinear Systems Including 
the Black Hole Equilibrium System, Comm. Math. Phys. 99 (1985), 563-591. 

[4] P. T. Chrusciel, "No Hair" Theorems — Folklore, Conjectures, Results, to 
appear in Proc. of AMS/CMS Conf. on DifF. Geom. and Math. Phys., Vancouver, 
August 1993, Eds. J. Bccrn and K. Duggal, Contemporary Mathematics, AMS, 
1994. 

[5] P. T. Chrusciel and R. M. Wald, On the Topology of Black Holes, ESI 

preprint, Vienna, October 1994. 
[6] P. Eberlein and B. O'Neill, Visibility Manifolds, Pacific J. Math. 46 (1973), 

No. 1, 45- 109. 

[7] F. ,]. Ernst, New Formulation of the Gravitational Field Problem, Phys. Rev. 

Letters 167 (1968), 1175-1178. 
[8] E. Heintze and H. C. Im Hop, Geometry of Horospheres, J. Diff. Geom. 12 

(1977), 481-491. 

[9] Y. Li and G. Tian, Nonexistence of Axially Symmetric, Stationary Solu- 
tion of Einstein Vacuum Equation with Disconnected Symmetric Event Horizon, 
Manuscripta Math. 73 (1991), 83-89. 

[10] , Regularity of Harmonic Maps with Prescribed Asymptotic Behavior and 

Applications, Comm. Math. Phys. 149 (1992), No. 1, 1-30. 

[11] S. D. Majumdar, a Class of Exact Solutions of Einstein's Field Equations, Phy. 
Rev. 72 (1947), 390-398. 

[12] P. O. Mazur, Proof of Uniqueness of the Kerr-Newman Black Holes, J. Phys. 
A 15 (1982), 3173-3180. 

[13] G. D. Mostow, Strong Rigidity of Locally Symmetric Spaces, Annals of Mathe- 
matics Studies, No. 78, Princeton University Press, Princeton, 1973. 

[14] A. Papapetrou, A Static Solution of the Equations of the Gravitational Field for 
an Arbitrary Charge-Distribution, Proc. Roy. Irish Acad. A51 (1947), 191-204. 

[15] R. Penrose, Some Unsolved Problems in Classical General Relativity, in Seminar 
on Differential Geometry, edited by S. T. Yau, Annals of Mathematics Studies 
no. 102, Princeton University Press, Princeton 1982. 

[16] D. C. Robinson, Uniqueness of the Kerr Black Hole, Phys. Rev. Letters 34 
(1975), 905-906. 

[17] R. SCHOEN AND K. UllLENBECK, A Regularity Theory for Harmonic Maps, J. 
Diff. Geom. 17 (1982), 307-335. 

[18] , Boundary Regularity and the Dirichlet Problem for Harmonic Maps, J. 

Diff. Geom. 18 (1983), 253-268. 

[19] R. Schoen and S. T. Yau, Compact Group Actions and the Topology of Mani- 
folds with Non-Positive Curvature, Top. 18 (1979), 361-380. 

[20] G. Weinstein, On Rotating Black Holes in Equilibrium in General Relativity, 
Comm. Pure Appl. Math. 43 (1990), 903-948. 

[21] , The Stationary Axisymmetric Two-Body Problem in General Relativity, 

Comm. Pure Appl. Math. 45 (1992), 1183-1203. 

[22] , On the Force between Rotating Co-Axial Black Holes, Trans. Amer. 

Math. Soc. 343 (1994), No. 2, 899-906. 



32 



GILBERT WEINSTEIN 



[23] , On the Dirichlet Problem for Harmonic Maps with Prescribed Singular- 
ities, Duke Math. J. 77 (1995), No. 1, 135-165. 

Department of Mathematics, University of Alabama at Birmingham, Birm- 
ingham, Alabama 35205 

E-mail address: weirLstei@math.uab.edu 



